Phase Estimation with Non-Unitary Interferometers: Information as a Metric 
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Determining the phase in one arm of a quantum interferometer is discussed taking into account the 
three non-ideal aspects in real experiments: non-deterministic state preparation, non-unitary state 
evolution due to losses during state propagation, and imperfect state detection. A general expression 
is written for the probability of a measurement outcome taking into account these three non-ideal 
aspects. As an example of applying the formalism, the classical Fisher information and fidehty 
(Shannon mutual information between phase and measurements) are computed for few-photon Fock 
and NOON states input into a lossy Mach-Zehnder interferometer. These three non-ideal aspects lead 
to qualitative differences in phase estimation, such as a decrease in fidelity and Fisher information 
that depends on the true value of the phase. 



I. INTRODUCTION 

Optical interferometers [l| and matter wave inter- 
ferometers [3] have been of great interest because of 
their practical applications in metrology. Interferome- 
ters have been used to measure such diverse quantities 
as electric, magnetic, and gravitational fields, gravita- 
tional waves @tJ], and there are plans to use them to 
test the theory of general relativity Q. Classical optical 
interferometers @ have been routinely used for sensing 
rotation in gyroscopic applications based on the Sagnac 
effect and experiments with Sagnac interferome- 

ters have been done with sin gle-p hotons [l2j |. with Bose- 
Einstein condensates(BEC) |13l - [l5| . and schemes using 
entangled particles have been proposed that are capable 
of Heisenberg limited precision measurements that scale 
as where N is the number of particles [l^. 

On a more fundamental level, there is interest in inter- 
ferometers because they are a vehicle to study the limits 
of precision of quantum measurements |17H19l | . Perhaps 
the simplest generic measurement problem consists of de- 
termining the relative phase shift between two arms of an 
interferometer from measurements made at the output 
ports of the interferometer [20l - |25| . This phase shift may 
be related to a classical external field incident on a phase 
shifter in one arm of the interferometer, in which case the 
interferometer can be used as a sensor of the field [26| . 
The determination of the phase shift is a specific exam- 
ple of the more general problem of parameter estimation, 
whose goal is to determine one or more parameters from 
measurements [27l - [33 |. 

Recently, there have been experimental demonstra- 
tions using entangled states to estimate the phase shift 
in one arm of a Mach-Zehnder interferometer |2ll [35l - l37| . 
Even more recently, the effect of losses on phase deter- 
mination was studied experimentally (ssj . In real ex- 
periments, there are three non-ideal elements of the in- 
terferometer system: state preparation [sgI Isoj , photon 
losses in the interferometer [385 and non-ideal photon- 
number detection [2l] . i37i] . Phase estimation has been 
theoretically investigated by separately takin g in to ac- 
count non-deterministic state preparation [28l434| . pho- 



ton losses in the interferometer itse lf 124 |25|. l4Cl| - |45l | and 
photon-number counting efficiency sj 4(i|. 

In this work, I write down a formalism that simulta- 
neously takes into account, in a unified way, all three of 
the non-ideal elements in experiments: non-deterministic 
state preparation, propagation through a lossy interfer- 
ometer, and imperfect state detection. Non-deterministic 
state preparation must be described by a density matrix, 
rather than by a pure state, thereby allowing for the finite 
probability of creating states other than intended. When 
the optical state is created, it enters the interferometer, 
where propagation may be non-ideal because photon ab- 
sorption and scattering can occur. Finally, when the opti- 
cal state leaves the interferometer, it enters the detection 
system, which may also be non-ideal: the state registered 
by the detection system may not be the true state that 
entered the detection system. 

Much of the previous work was focused on determin- 
ing the optimum measurements for determining phase 
and hence the quantum Fisher information was of pri- 
mary interest, because it gives a bound on the vari- 
ance of the phase associated with the optimum measure- 
ment LSSriSSj, i4^i49|] . In contrast, in this work I look 
at the information gain from specific, simple, photon- 
number counting measurements that can be easily imple- 
mented in the laboratory, and hence the classical Fisher 
information is the quantity of interest because it depends 
on the particular measurement that is performed. 

In Section |TT1 I briefiy review the theory of phase de- 
termination based on parameter estimation (Fisher in- 
formation) and on fidelity (Shannon mutual information 
between measurements and phase). In Section IlIIl I in- 
troduce an example of a non-ideal interferometer system, 
where state evolution is non-unitary. I write a statistical 
expression for the probability of measurement outcomes 
that takes into account the three non-ideal components 
of the interferometer system described above. I use this 
probability in the classical Fisher information in Eq. ^ 
and in the fidelity in Eq. ^ to analyze the determination 
of phase in a non-ideal interferometer system. As simple 
examples of the formalism, in sub-sections of Section Hill 
I look at few photon examples of non-deterministic state 
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preparation, propagation through a non-unitary (lossy) 
interferometer, and imperfect state detection. Finally, in 
Section IIV| I make some concluding remarks. My goal 
is to look at examples of few-photon states that can be 
implemented experimentally, with the hope that the ex- 
amples and method described here can be helpful for an- 
alyzing real experiments. Furthermore, in this work, I 
restrict myself to the simple case of non- adaptive mea- 
surements [soj . where the measurement is fixed before 
phase estimation. 



II. THEORETICAL BACKGROUND 

The accuracy of estimating a (single) one-dimensional 
parameter, (j), is described in terms of the classical 
Cramer- Rao bound [Slj, which gives a lower bound on 
the variance (Sip)'^ of an unbiased estimator of the pa- 
rameter (h: 
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Fcii(t>;M) 



(1) 



where Fci{(j); M) is the classical Fisher information given 
by MM 



pmp) 



(2) 



The classical Fisher information is described in terms 
of the conditional probability distribution, P(^|(/), p), for 
measurement outcome, ^, which can take one or more 
continuous values, or, one or more discreet values. If ^ 
takes continuous values, the sum over ^ in Eq. ^ is an 
integral. For the case of quantum measurements, these 
probabilities are given by 



P(e|<^,p)-tr /)^flo(0 =tr p„fl^ (0 



(3) 



where the state is specified by the Schrodinger pic- 
ture density matrix, p^, and the measurements by 
the positive-operator valued measure (POVM) in the 
Schrodinger picture, no(C)- The POVM are set a of non- 
negative Hermitian operators, M = {n(^)}, represent- 
ing a given physical measurement and so the expectation 
value of each operator no(0 is non- negative, and satisfies 
X]no('C) — I J where / is the identity operator. Alterna- 

tively, in the Heisenberg picture, the state is given by the 
density matrix, po, and probabilities of measurements by 
POVM, n<)4(^). Operators in the Schrodinger and Heisen- 
berg pictures, Os and Onift)), respectively are related by 
0^(0) = ;7t(0)O5 [/(</)), and U{(t)) = exp(-i0/i), where 

h is the infinitesimal displacement operator for the pa- 
rameter 0, satisfying 



(4) 



The quantum Fisher information, Fglcf)) is obtained by 
maximizing the classical Fisher information, Fci{(t>', M), 
over all possible measurements M, at a given value of 
(j). Braunstein and Caves [3l|, [s^ have shown that an 
improved lower bound is possible for the variance, {S(j)) , 
in terms of the quantum Fisher information: 
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(5) 



where Fq{(j)) is independent of the measurement M. The 
quantum Fisher information is defined by 



Fq (0) = tr p^Al 



(6) 



where the Hermitian operator, A0, is the symmetric log- 
arithmic derivative (S.L.D.), defined implicitly by 



-I- p^A^ 



(7) 



The right side and left side of the inequality in Eq. ^ 
is sometimes called the quantum Cramer-Rao bound, see 
also the work by Helstrom [2^ [2^ and Holevo [30j| , and 
discussion by Barndorff- Nielsen et al. [13, [13] ■ The ex- 
pression in Eq. ^ provides a bound on the variance of an 
unbiased estimator for an optimum measurement. How- 
ever, the theory does not give a procedure for determin- 
ing the optimum measurement. For one-dimensional pa- 
rameter estimation and for simple (non-adaptive) mea- 
surements, Barndorff- Nielsen and Gill have shown that in 
general the optimum measurement M will depend on the 
parameter (j), which is unknown prior to estimation [33| . 
Consequently, Barndorff-Nielsen and Gill have proposed 
a two-stage adaptive measurement procedure that will 
give 



Fei((/.;M) = FQ(0) 



(8) 



for optimum measurement M for all (fi. 

For the case of a pure state, \ipo), where the density 
matrix is po = IV'o) {'(pol, and where the path is generated 
by a unitary transformation, U{4>), the quantum Fisher 
information, Fq{4>), does not depend on (j) [H, [H, l49| . 
and is given by the fluctuations of the generator h by 
Fq {4>) = 4(A/i)^. Furthermore, Hofman has shown that 
for pure states having a path symmetry [s^l, the quan- 
tum Cramer- Rao bound in Eq. ^ can be achieved at any 
value of phase 0. The condition for optimal measure- 
ments for the case of pure states has also been investi- 
gated P]. 

In Section IIIII below, I consider the case of a lossy 
Mach-Zehnder interferometer, where the state evolution 
is effectively non-unitary, thereby leading to a classical 
Fisher information that depends on the true value of the 
phase <p. 

The above discussion of parameter estimation is based 
on classical and quantum Fisher informations, which are 
local descriptions of phase estimation, because they de- 
pend on the true value of 4>. Complementary to the above 
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local descriptions, is a global description given by the fi- 
delity [26]: 



H{M) = Y, I dcf,P{^\^,p)p{cl>) X 

C J — TT 



l0g2 



(9) 



where P(^|(/),p) is given in Eq. The fidelity, H{M), 
is the Shannon mutual information [sil . [s^ j between the 
measurement M and the unknown parameter (p. The fi- 
delity, H{M), gives the average amount of information 
(in bits) about the parameter (j) that can be obtained 
from the measurement M for one use (one measurement 
cycle) of the interferometer. The fidelity does not depend 
on (f) because it is an average over all possible phases 
and over all probabilities of measurement outcomes for a 
given POVM M. (For an alternative discussion of local 
versus global phase estimation, see Ref.[23|.) The fidelity 
also depends on the prior information about the param- 
eter (p through the prior probability distribution p((/)). 
Consequently, the fidelity characterizes the quality of the 
interferometer system as a whole, in terms of mutual in- 
formation between the measurement M and the param- 
eter (f). Note that the fidelity depends on the input state 
density matrix, p, and the measurement M, and there- 
fore can be used to optimize the system with respect to 
the input state and measurement. The fidelity is a mea- 
sure of the information that flows from the phase (j) to 
the measurements, which is analogous to a communica- 
tion problem where Alice sends messages to Bob. In the 
case of the measurement problem, quantum fluctuations 
in the initial state, in the channel (interferometer), and 
the type of measurement, determine the amount of infor- 
mation that is obtained about the parameter (j) from the 
measurements. The fidelity has been applied to compare 
the use of Fock states and NOON states when no prior in- 
formation is present about the phase [26] and when there 
is significant prior information about the phase (ssj . 

The complimentary measures of fidelity and Fisher in- 
formation may be contrasted as follows. Assume that I 
want to shop to purchase the best measurement device 
to determine the unknown parameter (jj. If I do not know 
the true value of the parameter (p, I would compare the 
overall performance specifications of several devices and 
I would purchase the device with the best over-all spec- 
ifications for measuring (j). The fidelity, H{M), is the 
over-all specification for the quality of the device, so I 
would purchase the device with the largest fidelity. After 
I have purchased the device, I want to use it to deter- 
mine a specific value of the parameter (f) based on several 
measurements (data). This involves parameter estima- 
tion, which requires the use of Fisher information, and 
depends on the true value of the parameter (j). 

Historically, the variance, (Scj))'^, of the estimated pa- 
rameter (j> has been discussed in terms of the stan- 
dard quantum limit, ScpsL = and the Heisenberg 
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FIG. 1. (Color) Interferometer system shown with three com- 
ponents: state preparation, interferometer, and detection sys- 
tem. 



limitjl, [H, [H, [13, 5(t)HL = where N is the num- 

ber of particles or quanta that enter the interferometer 
during each measurement cycle. The value 5(j) is presum- 
ably the width of some probability distribution, p(0|^, p), 
such as the distribution calculated from Bayes' rule, see 
Eq. (|25|) . Detailed calculation of p((/)|^, p) for a number of 
input states shows that these distributions have multiple 
peaks J6]. Consequently, rather than using the widths 
of these distributions as a metric for determining I 
use the information measures, Fisher information and fi- 
delity, which naturally handle distributions with multiple 
peaks. 



III. NON-IDEAL OPTICAL SYSTEM 

As described in the introduction, an interferometer sys- 
tem can be divided into three parts: state creation, state 
evolution through the optical interferometer, and detec- 
tion of the output state. In a real experiment, each of 
these three parts can be non-ideal, see Fig. [T] For exam- 
ple, I may want to create a quantum state [■i/;™) as input 
into the interferometer. However, instead, the resulting 
state may be a mixture of states, each with some prob- 
ability, Ps('/'fe")j for fc = 1, 2, • • • . Such a quantum state 
is described by the density matrix p: 



E 

k 



(10) 



The state p is then input into the interferometer, where 
there may be absorption and scattering of photons. For 
example, a two-photon state may enter the interferometer 
and a one-photon state may exit the interferometer, be- 
cause one photon was absorbed inside the interferometer. 
Alternatively, a two-photon state may enter the interfer- 
ometer and a three-photon state may exit the interferom- 
eter, due to light scattering into the interferometer from 
the environment. I can describe these processes gener- 



ally by a transfer matrix, Pj {ip" Itpl. 



I, which gives 
*) to exit the 



the conditional probability for state 
interferometer given that state IV'™) entered the interfer- 
ometer. The transfer matrix, P/(i/'°"*|'0™, (/)), is general 
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enough to describe non-unitary propagation of the quan- 
tum state through the interferometer, and so can take 
into account losses and scattering. Note that the trans- 
fer matrix may depend on the state of the interferome- 
ter, which I specify here by single parameter 0. Finally, 
the detection of the quantum state that leaves the inter- 
ferometer can be non-ideal. For example, the detection 
system may register a measurement ^, when state IV'i'"*) 
enters the detection system, whereas the true state that 
entered the detection system was I can represent 

such an imperfect detection system by the conditional 
probability PcCCIV'j"*) which gives the probability for 
making a measurement ^ when state entered the 

detection system. Note that in general this probability 
may or may not depend on (/), a parameter describing the 
state of the interferometer. For a non-ideal interferome- 
ter system, t he p robability of obtaining a measurement 
f is given by [5^ 

j k 

(11) 

where we must have each of the three probabilities sum 
to unity: 

E ^^^(^r) = 1 (12) 

k 

E^/('^ri^r,0) = i (13) 
Y,PDmT'A) = ^ (14) 

i 

Equation pT|) is a general statistical relation for the 
probability of obtaining a measurement outcome ^ for 
given phase shift 0, taking into account the three non- 
ideal aspects of interferometer systems. Note that 
Eq. pT|) is sufficiently general that it can be applied to 
the case where states are represented by density matrices. 
In this case, in Eq. pT|) we can make the replacements 
■0™ -)■ and tA""* /Oj"*, where and are a 
set of input and output density matrices labeled by inte- 
gers j, fc = 1, 2, • • • . In order to compute the probabilities 
of measurement outcomes, P{^\4>), Eq. ([TT|) must be aug- 
mented by a detailed model of input and output states. I 
give several examples of applying Eq. (jlip in the sections 
that follow. 

The probability of measurement outcome, given by 
Eq. (|lip. enters into the Fisher information and into the 
Shannon mutual information, in Eq. ([2]) and Eq. ([9]), re- 
spectively. In the next three subsections. A, B, and C, 
I give examples of the effects of non-deterministic state 
preparation, state evolution through an interferometer 
when absorption is present, and imperfect output state 
detection, respectively, using Fisher and Shannon mutual 
informations as metrics of performance of the interferom- 
eter. 



A. Non-Deterministic State Preparation 

Consider an optical interferometer system that has 
non-deterministic state preparation, but has no losses 
in the interferometer and has perfect state detection. 
When I try to prepare a certain quantum state for in- 
put into the interferometer, there is always a non-zero 
probability that another state than intended will be pre- 
pared. This non-deterministic state preparation is ex- 
pressed by a density matrix for the input state, which 
assigns probabilities for creating various quantum states, 
see Eq. pHI) . Since state detection is assumed perfect, 
Pd(^|-0°"*, (/i) — 1 when the measurement f corresponds 
to the true state that entered the detection system, ?/;°"*, 
and otherwise ^^(flV''"'*, 0) = 0. 

A general interferometer with no losses is characterized 
by a unitary scattering matrix, Sij{4>), that connects the 
Np input-mode field operators, di, to the Np output field 
operators, Pi. 

A = E ^^j ('^) = UHcj))a,U{(b) (15) 

where U{(j)) is a unitary evolution operator, i,j — 
1, 2, • • • , Np, and (f) is one or more parameters (e.g., phase 
shift) that describe the state of the interferometer. 

For simplicity, I consider a Mach-Zchnder interferom- 
eter, with no losses, with input ports labeled, "a" and 
"b", and output ports, "c" and "d", having a scattering 
matrix 

Q 1 f -lil + e'-^) (-l + e'"^) \ 

S^A^) = ^[ + z(l + e^^) j ^16) 

where cti = {a, 6} and /3i — {c, d}. 

The probabilities, P/(V'°"*|V'r' </>)' t^iat relate the in- 
put state ■f/'fc" to the output state -0°"* of the interfer- 
ometer are given in terms of the projection operators 
{nc,nd)- 

p/(V'-*i0r,'/') = (V'rl inc,n,) |0r) (17) 

where the output state ■0^"* is specified by two inte- 
gers, {ric, rid}, giving the photon numbers output in ports 
"c" and "d". In terms of the unitary evolution opera- 
tor, U{<p), the output state in the Schrodinger picture is 
|.0°"*(0)) = [/(0)|0™), where jV'™) is the Heisenberg 
picture input state. In Eq. (jlip . the sum over k is now 
a double sum over all non-negative values of the two in- 
tegers ric and n^. For a generic Mach-Zehnder interfer- 
ometer, with input ports "a" and "b" , and output ports 
"c" and "d" , the projective operators are ^ 

n^(n„n,) = ^^ (ct)"- (dt)". |0) (0| {cr (dT" 

(18) 

where the vacuum state |0) = |0)^ (g) |0)j,. 
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FIG. 2. (Color) Lossy Mach-Zehnder interferometer is shown, 
with input modes ai, a2, vi, V2, and output modes bi, &2, di 
and d2- Here modes vi and V2 have vacuum input and the 
output modes di and d2 are loss channels in each arm. 



As an example, consider the simplest case of input 
given by a mixed state represented by the density ma- 
trix 



p = Po |0) (0|+Pi 1 10) (10| 



(19) 



where Pi is the probability for 1-photon input into port 
"a" and vacuum input into port "b" , and Pq is the prob- 
ability of vacuum input into both ports "a" and "b" , and 
Po Pi = 1. In Eqs. (Uni) and (HI]), Ps{^'") = Po when 
IV'™) = |0) and PsiV) = Pi when jV-*") = |10). 

Equation p?|) for the interferometer transfer matrix 
can then be written as 



P{nc,nd\(j),p) = tr (^pfl^, {nc,nd) 



(20) 



where the input state is represented by the density ma- 
trix p. It seems that there can be only two possi- 
ble measurement outcomes, ^ = {ncn^} = {1,0} and 
= {nc,nd} — {0,1}. However, the probabilities for 
the two measurement outcomes do not sum to unity be- 
cause, P(1O|0, p) + P{01\(p,p) = Pi. Therefore, there is 
a non-zero probability of an inconclusive measurement 
outcome associated with the probability Pq of vacuum 
injected into both input ports "a" and "b" . I introduce 
an inconclusive measurement operator, n0(i), so that 
the sum of the three operators is equal to the identity 
operator /: 



Il,p{lO)+U^ (01) 



(21) 



Using Eqs.([T8l)-(l2T]), the probabilities given by Eq. (ITT1) . 

P(^|(/i), for measurement outcomes ^ are given by 



(22) 



P(1O|0) = Picos^ ( ^ 

P{Ol\(j)) = Pisin^ 
Pii\q}) = 1-Pi 



(23) 
(24) 

where P{i\(f)) is the probability for an inconclusive mea- 
surement outcome. 



Using Bayes' rule and Eq. (|TT|) . we can write the con- 
ditional probability distributions for the phase, p((/>|C), 
given measurement outcome, ^, as 



Pirn - 



p(ci0)p(0) 



(25) 



/ Pimpmd(^' 



1 2 1 
- COS 

TT 


it) 


(26) 


1.2/ 

— sm 

TT \ 


;i) 


(27) 


1 

2^ 




(28) 



where p{(j)) is the prior probability distribution specifying 
our prior information about the phase (f). Assuming no 
prior information about the phase, p('/>) = l/(27r), and 
using Bayes' rule in Eq. ([25]), the conditional probabil- 
ity distributions for the phase for a given measurement 
outcome are given by 



Pirn) 

p(0|Ol) 
p{(p\i) 



where f — (ric, rid)- As we would expect, for an inconclu- 
sive measurement outcome the phase probability distri- 
bution, p(0|i), is flat since an inconclusive measurement 
result cannot be used to distinguish different values of 
the phase. Note that the probabilities for the three mea- 
surement outcomes in Eq. sum to unity. When 
Pq = 0, or equivalently Pi = 1, state preparation is deter- 
ministic, and the probability for an inconclusive measure- 
ment outcome is zero. In this case, the probabilities for 
measurement outcomes (10) and (01) reduce to the values 
for the case of a single-photon input state created with 
probability unity. Note that this single photon input case 
corresponds to the case of a classical interferometer fed 
by a laser in one input port. Larger photon-number in- 
put states have measurement outcome probabilities that 
differ from the probabilities for a classical interferometer. 

For the input state in Eq. ([T51) . the classical Fisher 
information, defined in Eq. ([2]), is given by 



F{cp) = Pi 



(29) 



where I dropped the subscript cl on the classical Fisher 
information, a convention that I follow in rest of this 
work. According to the Cramer- Rao bound in Eq. ([l}, 
when the probability of creating a single photon ap- 
proaches zero. Pi — >■ 0, the variance becomes ar- 
bitrarily large, because the probability for inconclusive 
measurement outcomes approaches unity. 

Assuming no prior information about the phase, there- 
fore taking p{(f) = l/(27r), the fidelity (Shannon mutual 
information) of the system defined in Eq. ^ is given by 



iI(M).Pi(^-l 



(30) 



Similar to the Fisher information, the fidelity H{M) also 
approaches zero when the probability Pi of having one 
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photon in the input of each shot approaches zero. The fi- 
dehty is the amount of information (in bits) that is gained 
on average about the phase from a single use of the in- 
terferometer, averaged over all possible phase values (j). 



B. Lossy Mach-Zehnder Interferometer 

Next, I consider an interferometer with absorption 
losses — so state evolution is non-unitary. I assume that 
state preparation is deterministic (ideal) and that state 
detection is perfect (no errors). Equation (fTTj) is gen- 
eral enough to describe processes other than losses in the 
interferometer, such as photons scattering into the inter- 
ferometer from the environment, in which case there are 
more photons leaving the output ports than entering the 
input ports. However, in what follows, I restrict myself 
to simple absorption in the interferometer. I model losses 
in each arm of a Mach-Zchnder interferometer by insert- 
ing two beam splitters, S3 and S^, one in each path, see 
Fig.[2j While a lossless Mach-Zehnder interferometer has 
two input and two output ports, a general lossy Mach- 
Zchnder interferometer can be represented by four input 



and four output ports, see Fig.[2l I label the input modes 
as ai, a2, wi, and W2, where vi, and V2 have vacuum input 
and I label the output modes as 61, 62, and d2, where 
di , and ^2 are the modes where probability amplitude is 
"dissipated". I take the phase shifts at the two mirrors, 
All and M2 to be equal to tt. Furthermore, I assume 
that the interferometer is balanced, so that path lengths 
satisfy. 



L = li+l:i + k = l2 + U + k 



(31) 



see Fig. [21 A calculation gives the input and output 
modes related by the 4x4 unitary scattering matrix 




Sij (0) 




(32) 



where the phase-dependent scattering matrix is given by 



2^ 



J4> 



2*^ 



1 — r'^ 

y 



7^' 



«(0+%^) 



(v^ 



1 _ ^2 

y 



%/2 ^ 



^/2^^'^ 







(33) 



It is easy to check that the scattering matrix is unitary, 
S = I where / is the 4x4 unit matrix. The parame- 
ters, and Vy, are the reflection amplitudes for beams 
splitters 5*3 and ^4, respectively, and they represent the 
strength of the loss or dissipation, see Fig. [5J When the 
system is considered in terms of two input modes, ai and 
a2, and two output modes, bi and 62, the evolution of the 
input state in not unitary. 

The 4x4 unitary scattering matrix Sij{(t>) has some 
simple properties. The case when — ry — corre- 
sponds to no dissipation. In this case, the 4x4 S- matrix 
reduces to two diagonal 2x2 blocks. The upper left 2x2 
block couple modes ai and 02 to modes bi and 62, and 
this 2x2 block (up to a phase) is given by Eq. (|16l) , which 
is the scattering matrix for the Mach-Zehnder interferom- 
eter with no losses. For this case of no loss, in Eq. ([55)1 
the lower right 2x2 block couples the dissipative modes, 
di, and (i2, to the vacuum modes, vi, and V2- 

The case — ry — 1 corresponds to maximum dissi- 
pation, and the 4x4 S- matrix again decouples, into two 
off-diagonal 2x2 blocks. The upper right 2x2 block cou- 
ples the two vacuum modes, vi and V2, to the two out- 
put modes, 61 and 62- The lower left 2x2 block of this 



S-matrix couples the loss modes, di, and c?2, to the input 
modes ai and 02. For this case of maximum dissipation, 
the input modes, ai and 02, are decoupled from the out- 
put modes, hi and 62. 

The probabilities for various measurement outcomes 
^ = (n, m) are given by the analog of Eq. ([3]): 

P{n,m\(j),p) ^tT(^Pofl^{n,m)j (34) 

where n and m are the number of photons leaving ports 
bi and 62, respectively. The trace is over the com- 
plete space of four direct product Fock basis states, 
hi)ai l"-2)a2 ® l"3)t,i ^4)^^- Thc input state den- 
sity matrix, po, is defined in terms of a sum of products 
of creation operators d\, d^, v\ and acting on the vac- 
uum |0) = |0)^^ (g) |0)^^ (g) |0)^^ (g) |0)^^ and has the form 

po — Cnm \nmOO) {nmOO\ (35) 

where I use the short-hand notation 
|nmOO) = |n)„^ (g |m)„^ (g |0)^^ (g |0)„^. I am using 
the Heisenberg picture, so the input density matrix 
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Po is independent of time (phase), while the operators 
(n, m) evolve in time (phase) and so they depend on 



The projective measurement operators are given by 



(n, m) 



nl m 



k,l=Q 



(^0" (^0'" i^^'T (^"0' i'^^ (^0" i^'T i'^'T 



(36) 



where b\, feji '^i: ^nd are creation operators for the 
output modes bi, 62, di, and ^2, respectively. The sums 
over k and / take into account the probabilities for losing 
photons into ports di and o?2- I use the short-hand nota- 
tion for the input modes {ai} = {ai,a2,vi,V2} and for 
the output modes {Pi} = {61, 62, lii, ^2}, see Eq. (ITS|) . 

Note that, even though dissipation is being modeled, 
it is easy to check that there is global photon number 
conservation, 



(37) 



1=0 



1=0 



since the S-matrix in Eq. ([33]) is unitary. The measure- 
ment outcomes, ^ — (n, m), are specified by two integers, 
which label the number of photons that are output at 
ports bi and 62, see Fig. [2] So the probability of output 
state, as given in Eg. (13^ . is expressed by the two 

integers n and m, see also Eq. ((TT)) . Note that in general 
the number of photons n -I- m in the output state, V'J"*, 
is not equal to the number of photons in the input state, 
ipin, because 



a[ai + ala2 7^ b[bi + 6^62 



(38) 



However, the sum of probabilities for all possible mea- 
surement outcomes is unity: 



^ P{n,m\(j),p) = 1 



(39) 



m,n=0 



where P{n,m\(j), p) is given by Eq. p4l) . and is a result 
of 



^ (n, m) = / 



(40) 



For the case of a pure input state, \ipin), it is convenient 
to define the operators 



N (n, TO, fc, /) 



1 



n! to! k\ II 



(h) (52) (J.) (4) 



and the probabilities in Eq. ([3^ are then given by 



(41) 



P(n, to|(/), ^„) = ^ I (0| N (n, m, k, I) | V'™) (42) 



fe,;=i 



r 



Equation (p4)) . or Eq. (1421) for the case of pure states, 
defines a unitary mapping, |nmOO) — ?• jn'TO'fc^), be- 
tween interferometer input states, jnTOOO), and output 
states, \n'm'kl), because the photon number is conserved: 
n + m — n' + ra' + k + l, see Eq. ((37)) . 

If we restrict our attention to measurement out- 
comes projected onto the Hilbert subspace with basis 
W)b ® \m')^ , Eq. dSll) or Eq. ^ defines the non- 



unitary mapping £: 



® m 



(43) 



since photon number is not conserved: n-\-m ^ n' -I- m', 
see Eq. (p8)) . which represents losses in the Mach-Zehnder 
interferometer. Here the Fock states |?^')f)^^ * = li2, 
are created from the vacuum state, |0)b. , by application 
of creation operators &| in the usual way. The mapping 
E depends on two parameters, and rj,, which specify 
the strength of the dissipation or losses in each arm of 
the interferometer. In the limit of no dissipation, when 
r,j. — ^ and Vy — 0, the mapping E becomes a unitary 
transformation and n -V m — ri! ~\- vn! . 

In what follows, I use the short-hand notation \nm) for 
the input state InmOO) = |n)ai ® |™)a2 ® |0)t)i ® \^V2- 

In the next two subsections, IIII CI and IIIIDl I discuss 
the Fisher information and the fidelity (Shannon mutual 
information) for specific cases of few-photon Fock state 
and NOON state input into the lossy Mach-Zchndcr (MZ) 
interferometer . 



C. Fock State Input into Lossy MZ Interferometer 



Consider the A^-photon Fock state 



1 / 

IV^at) - -j= (at) |0) = lA^ooO) s |iVO) (44) 



input into a lossy Mach-Zehnder interferometer given by 
the scattering matrix in Eq. (|33)) . The probabilities for 
measurement outcomes ^ = (n, to) are given by: 



P(n,TO|(/),V'Ar) 



iV! 



N N 



EE 



n\m\ ^-^ ^-^ kill 

k=o 1=0 



on Qm Qk QL ^ 
'-'ll'-'21'-'31'-'41 "n 



+m+k+l,N 



(45) 



where Sij are the matrix elements of Eq. (|33]) and (5„j_„ 
is the Kronecker delta function. 
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A direct calculation of the classical Fisher information 
for the A^-photon Fock state input, F/v (</>), gives 



(46) 



where Fi{(j)) is the classical Fisher information for one- 
photon input, given in Eq. ()53p . This shows that for 
the lossy MZ interferometer with A^-photon Fock state 
input, the standard deviation (Scj)) scales as 1/a//V- From 
another point of view, since the Fisher information is 
additive for independent events, the A^-photon Fock state 
acts like N independent 1-photon states. When — ry, 
then Fi{(f)) — 1, and the A^-photon Fisher information 
becomes Fj^^cj)) = N. This means that dissipation in the 
(non- unitary) interferometer has the effect of introducing 
a phase dependence into the Fisher information, see the 
discussion below. 



1. 1-Photon Fock State Input into Lossy MZ Interferometer 

As the simplest example of the effect of dissipa- 
tion, I consider the 1-photon Fock state input into the 
lossy Mach-Zehnder interferometer with scattering ma- 
trix given by Eq. (155)) 



\^rn)=a\\{)) = 11000)= 1 10) 



(47) 



where I use the short-hand notation |10) for the input 
state 1 1000). The probabilities for measurement out- 
comes are given by P(n, m|0, ?Ai„), where n,m specify 
the photon numbers output in port hi and 62, respec- 
tively, see Fig. m The probabilities P{n,m\(t),iJ;in) for 
the three measurement outcomes are: 



1^2 



Fy) COS(f> 



(2^rl-rl + 2y'(l-r-i)(l-r-2) cos 0) 



P{W\(I>, 10) 

P(01|<^, 10) 
P(00|<?!>, 10) 

(48) 

The probability P{00\(j), 10) is associated with an incon- 
clusive measurement outcome, since for this case zero 
photons leave the output ports, i.e., the photon that en- 
tered in port "a" was absorbed in the interferometer, or 
more precisely the photon was output in either port di 
or ^2- 

From Bayes' rule in Eq. ([25]). the phase probability 
distributions, p{(l>\mn,ipin), for input state IV'm) given 
in Eq. (H71) . are given by 



P(0|1O,1O) = 27 2-r^-r^ 



p(0|Ol,lO) 
P(0|OO,1O) 



^ 2-rl-rl+2^il-r{){l-rl) cos,/, (49) 

J_ 

2n 



When ^ ry, there is a loss of contrast in the phase 
probability distributions p{4>\mn, 4'in), see Fig. [3] When 
the absorption probabilities are the same in both arms, 
in the limit = r^, the phase probability distributions 
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FIG. 3. (Color) For the 1-photon input state |10), given 
by Eq. (|47p . the probability distribution for the phase, 
p((;/)|10, 10), in Eq. (|49[) is plotted for absorption = and 
=0.0, 0.90, 0.95, 0.98, and 1.0. As -^1.0, the probabil- 
ity distribution becomes flat and does not distinguish between 
difl'erent phase values. 



in Eq. (j49| reduce to the case of a single photon input 
without losses in the interferometer, which are given in 
Eq. (|26l) - ([28)) . with trivial phase change given by the re- 
placements sin — )> cos. The phase probability density, 
p((/)|00, 10), is associated with the inconclusive outcome, 
p((/)|i). It is a remarkable feature that for equal loss 
in both arms, r^, = Vy, the phase probability densities, 
p((/)|mn,-0i„), in Eq. (j49|) do not depend on the size of 
the loss, Tx- However, there is loss of information with 
increasing absorption, r-x and Vy, which is reflected in the 
information measures, see below. 

The effect of equal dissipation in both arms of the in- 
terferometer is the same as the effect of non-deterministic 
state preparation, specified by input state characterized 
by a density matrix in Eq. ()19p . However, when the dis- 
sipation in both arms is not equal, say for Tj, = and r^, 
is finite, the phase probability distributions show a loss 
of contrast, see Fig. [3] This feature may be useful in 
applications to null-type measurements. 

In the discussion that follows, I assume no prior infor- 
mation on the phase, so I take p((/)) = l/(27r). When 
there is no loss in the interferometer, r^; = Tj, = 0, 
the Shannon mutual information (fidelity) as defined in 
Eq. © is a constant: 



(50) 



When the losses in both arms are equal, = Ty, we have 
the exact result 



(51) 



For general values of and r^, the expression for 
H{M) is large and complicated, but for small Tj, <C 1 
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FIG. 4. (Color) For 1-photon Fock input state, |10), given by 
Eq. (|47ll . the fidelity (Shannon mutual information) H{M) 



seems linear in 



for 



however, this is not true for 



general values of and ry , see also Fig. (5] 




FIG. 5. (Color) For 1-photon Fock input state |1000) = |10), 
given by Eq. (|47|) . the fidelity (Shannon mutual information) 
H{M) is plotted as a function of loss parameters and r-y. 



and Tj, ^ 1, I can expand it in a power series, 



For the 1-photon Fock state (given in Eq. (|47])) in- 
put into the lossy Mach-Zehnder interferometer given 
in Eq. (p3| . the classical Fisher information (defined by 
Eq. ^) is given by: 



2 1- 



rl) sin^ 



(2-r2-r2)^-4(l-r2)(l-r2), 



" " " 32(</,) 

(53) 

see plots in Fig. |6l From these plots, it is clear that for 
a lossy interferometer, the Fisher information depends 
strongly on the true value of (j). Through the Cramer- 
Rao bound in Eq. (IT|), this translates to a dependence of 
the variance ((50)^ on the true value of 0. 

Therefore, Fisher information Fi {(/)) for 1-photon Fock 
state input into a lossy interferometer is qualitatively dif- 
ferent than for an ideal interferometer without loss, see 
Eq. (15^ for the case = Vy = 0. As a consequence of 
Eq. (P5)) . the Fisher information for A'^-photon Fock state 
input into a lossy interferometer is qualitatively different 
than for a lossless interferometer, where it is a constant 
given by Fn{(P) — N. Specifically, for a lossy interferom- 
eter the Fisher information depends on the value of the 
true phase, see plots in Fig. |6l For values of the phase 
given by = and cj) = ±7r, the Fisher information van- 
ishes for Fock state input, independent of the value of 
the dissipation parameters, and ry, see comment [59[. 
According to the Cramer- Rao bound, the variance (i50)^, 
is large for values of true phase near = and (f) — ±7r. 
However, when there is no dissipation, — Vy — 0, 
the Fisher information is independent of (j) and so is the 
bound on the variance (i5<^)^. The presence of dissipa- 
tion in the interferometer introduces a dependence of the 
variance, ((5(/>)^, on true phase 4>. The exception to this 
is when 



where Fi (cj)) reduces to Fi (0) ~ 1 



and is independent of ' 



H{M) 



1 

hi2 



1 / 2 , 2\ 

2 V'^ + ^y) 



+ 0{rt) + 0{rt) (52) 



where I dropped fourth order terms in and ry. When 
r^ — ry, from Eq. (f52|) . we may expect the fidelity (Shan- 
non mutual information) H{M) to be quadratic in r^, 
however, this is not true for general values rx and ry, see 
Fig. m which shows H{M) vs. for the case r^ = ry 
and for ry — 0. 

In Figure [5] the fidelity (Shannon mutual information) 
H{M) is plotted as a function of the dissipation param- 
eters, rx and ry. When the dissipation in either arm is a 
maximum, r^. = 1 or = 1, the fidelity H{AI) — 0, indi- 
cating that we obtain zero information from each photon. 



2. 2-Photon Fock State Input into Lossy MZ Interferometer 

Consider now the 2-photon Fock state input into the 
lossy Mach-Zehnder interferometer with S-matrix given 
by Eq. 



-^(al)'|0) = 12000)^120) (54) 



IV'm) 



where I use the short-hand notation |20) for the state 
1 2000). The probabilities P(m, n\(j), ipin) for the six mea- 
surement outcomes are given by 
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FIG. 6. (Color) The Fisher information F{(j)) for 1-photon input state |1000) (short-hand notation |10)), is plotted as a function 
of Tx and ry, for different values of = 0, 0.125, n/2, 3.0, 3.1, vr, left to right in top row and bottom row. 



P(20|(;/>,20) = 


1 

16 




6 — 6r^ 


~&rl + A^{r%~l) (r2-l) (r^ + - 2) 


cos(0) + 2 (r^ 


~ 1) {"^l ~ 1) cos(20) + Arlrl + rl + ry 


P(O2|0,2O) = 


1 

16 




6 — 6rl 


-6r^-4y'(r2-l) (r2 - l) [rl + rl - 2) 


cos(0) + 2 (r^ 


- l) (ry - l) cos(20) + Arlrl + rt + r^ 


P(11|0,2O) = 


\ 


2 - 2r^ - 


2rl + r* + - 2 (1 - rl) (l - r^) cos (20)] 




P(1O|0,2O) = 




rl + r-y) 


2-r^-r^-2^(l-r2)(l-r2) cos (0) 






P(O1|0,2O) = 




7-^ + rl) 


2 - - + 2y'(l - rg) (l - rg) cos (0) 






P(OO|0,2O) = 




r^, + ry) 









The sum of the probabilities for the six possible mea- 
surement outcomes in Eq. (j55p is unity. Since a two- 
photon Fock state has been used as input, the probabil- 
ity that no photon was absorbed is equal to the the sum 
P(2O|0) + P(O2|0) + P(11|0) = i (2 - r2 - rl)\ while 



the probability that exactly one photon was absorbed is 
P(1O|0) + P(O1|0) = i {2rl + 2rl - 2rlrl - - r^). 

The conditional probability distributions for the phase, 
p((f>\mn,'ipin) , with input state \4>in) in Eq. ((54)) and 
measurement outcome ^ = (m, n) are given by 
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p(0|2O,2O) 

p(</.|ll,20) 
p(0|lO,2O) 

p(0|Ol,2O) 
p(</.|00,20) 



1 

2^ 
1 

2^ 
1 

2^ 
1 

2^ 
1 

2^ 
1 

2^ 



4^(l-r2)(l-r2) {rl + - 2) cos((/)) + 2 (r^ - l) (r^ - l) cos(20)" 
^ 6-6(r2+r2)+r4+r4+4r2r2 



4^(1 - rg) (1 ~ rg) (rg + rg - 2) cos((/)) + 2 (rg ~ l) (r.^ - l) cos(20) 
6-6(r2+r2) + r4+r4+4r2r2 

2 (1 - r2) (1 - r2) cos(2^ ' 



1 - 



1 + 



2-2(r2+r2 


) + ^^ 


I ^4 

T- / J, 


2y(l-r2)(l 




C0S((/)) 


2~r2 


— r2 






-i) 


cos(0) 




— r2 





(56) 



In the limit of high photon absorption, — 1 and fidehty is given by 



ry = 1, the phase probabihty densities p{4>\mn, 20) = 
l/(27r) for all measurement outcomes, so that different 
phase values are not distinguishable because of photon 
losses. 

A remarkable property of the phase probability distri- 
butions in Eq. ([56]) (similar to that of Fock state input in 
Eq. (Uni)) is that for equal losses in both arms of the in- 
terferometer, r,j; — Ty, the phase probability distributions 
are independent of the magnitude of the loss , having 
values 



where I have written the values of the functions (from 
top to bottom) in Eq. as components of a vector left 
to right in Eq. (|57l) . One may think that for — Vy 
the contrast in the photon number measurements is lost, 
however, setting = ry in Eq. and expanding for 
small <C 1 leads to 



/ P{20 
P{02 
P{11 
P(10 
P(01 

V p(oo 



10,20) \ 
10,20) 
10,20) 
10,20) 
10,20) 
10,20)/ 



sm 



£ - 2r2 sin* £ 



cos ^ 



9 • ? 

rt sm 



cos(0)) 
cos(0) + 1) 
V ^ 0(4) 

(58) 

which shows that the probabilities become vanish- 
ingly small for the measurement outcomes P(1O|0, 20), 
P(O1|0, 20), and P(OO|0, 20), which correspond to one or 
both photons being absorbed, but contrast for different 
measurement outcomes is not lost. 



For the case of equal loss in both arms. 



with 



no prior information on the phase, p(0) = l/(27r), the 



H{M) = 



41n2 



9 / 11 

4 In 2 - 3 In 3 -I- 2r^arctanh — 

43 



(59) 

For the case of small (but not equal) losses in both 
arms, I and ry <Si 1, the fidelity is given by 



H{M) 



8-41n2-31n3 



31n3 



1 



41n2 ln2 



2 2 

r r 

X y 



41n2 

l + ln2 — ln3\ ^ , .4 ^ , ^4 /^„n 
^ ) + O (r^r + O (ry)^ (60) 



where I have dropped terms of fourth order in rx and ry . 
Equation (jBOj) gives the dependence on the dissipation 
parameters, r^ and ry, of the information gain about 0, 
for single use of the interferometer, when there is no prior 
information about 0. 



D. NOON State Input into Lossy MZ 
Interferometer 



Next, I consider the iV-photon NOON state input into 
the lossy Mach-Zehnder interferometer with scattering 
matrix given by Eq. p3|) : 



1 



NOON) 



(«I)%(4) 



N 



|0) 



V2 



[|iV000) + I07V00)] 



(61) 
(62) 
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For the input state in Eq. (|62|) . the probabihty for mea- 
surement outcome ^ — (n, m) is given by 



P{n,fn\(j),'ipNooN) 



2nlml k\l\ 

k=0 1=0 



5' n cm nk qL 
I 11'^21'-'31'-'41 



on QTn nk qI | 
'Jl2'-'22'-'32'^42 



Sn 



+m+k+l,N 



(63) 



amount of dissipation, r^, for Fock states the fidehty in- 
creases with input photon number N. The fidehty for 
2-photon NOON state input is exactly zero for all values 
of dissipation because this state carries no information 
about the phase in a Mach-Zehnder interferometer, see 
the discussion below. 



where n and m are the number of photons output in 
ports bi and 62, respectively. Using the Fisher informa- 
tion, I compare how well Fock states and NOON states 
perform in the presence of absorption losses. In Fig. [3 
I plot the classical Fisher information for TV =3, 4, and 
5 photon Fock states and NOON states, plotted vs. r^. 
for the special case where — Ty. The plots show that, 
for equal dissipation in both arms, and for equal photon 
number, Fock states perform better for phase estimation 
than NOON states, for the same amount of dissipation r^, 
see Eq. ((T]). 

For NOON states, the Fisher information vanishes at 
the phase values: (/) — 0,±7r/2,±7r. While for Fock 
states, the Fisher information vanished only at </> = 0, ±7r. 
For close to these values, phase estimation may have 
large standard deviation, see Eq.([T]). 

Figure [5] shows the classical Fisher information for the 
case where Fock and NOON states are input into a Mach- 
Zehnder interferometer with small dissipation (losses) 
and when the losses are not equal in both arms. Gen- 
erally, Fock states perform better (have larger Fisher in- 
formation) for all values of dissipation Vx except at the 
very highest values of r^; ~ 0.95. The comparison is made 
at a true value oi (j) = 7r/4, where the Fisher informations 
do not vanish. 

Figure [3] shows a plot of the classical Fisher informa- 
tion for Fock and for NOON states for = 3, 4 and 5 
photons for the case of large dissipation in one arm of 
the Mach-Zehnder interferometer, Ty = 0.9. The com- 
parison is complicated, since for the N —i photon case, 
Fock states perform better than NOON states for large 
dissipation ^ 0.8, whereas the situation is reversed for 
small dissipation ^ 0.05. 

In Figures [S] and ini the comparisons are made at a true 
value of phase = 7r/4, where the classical Fisher infor- 
mations (for Fock and NOON states) do not vanish. When 
dissipation is present, the classical Fisher information has 
a complicated behavior as a function of the true phase 
4>, see Fig. [TU] This shows that phase estimation using 
simple photon counting is a sensitive procedure whose 
accuracy depends on the true value of phase. 

The classical Fisher information for a lossy Mach- 
Zehnder interferometer depends on the true value of the 
phase 4>. The fidelity (Shannon mutual information) is 
an information measure that averages over all phases, for 
prior information given by p{4>), see Eq. Figure [TT] 
shows a comparison of the fidelity versus dissipation 
for Fock and NOON states for equal dissipation in both 
arms, = Ty. The fidelity of 1-photon Fock and NOON 
states is equal, see the discussion below. For a given 



1. 1-Photon NOON State Input into Lossy MZ 
Interferometer 

Consider now the 1-photon entangled NOON state: 



v2 



1 
71 



[|10) + |01)] (64) 



where again, I use the short-hand notation |10) for |1000) 
and |01) for |0100). The probabilities for the measure- 
ment outcomes for this input state are given by Eq. (j48p 
with the replacement cos (p — >■ sin (j). Similarly, the phase 
probability distributions, assuming no prior information, 
p(0) ~ l/(27r), are given by Eq. (l49l) with the replace- 
ment cos(j) — )■ sin0. The fidelity for this input state 
is the same as for the 1-photon Fock state, given by 
Eq. ((50)) - ([52)1 . Therefore, according to Shannon mutual 
information (fidelity), the presence of entanglement in 
the 1-photon NOON state has not improved the informa- 
tion on the phase. 

The Fisher information for this entangled state is given 
by the 1-photon Fock state Fisher information in Eq. ([55)) 
with the replacements 0—5- ^ ~ cf). Therefore, the en- 
tanglement simply has the effect of changing the phase 
of the classical Fisher information. This phase change 
changes the places where F(0) = 0, which, for this en- 
tangled state, is now (p = ±7r/2. Comparison of the 1- 
photon Fock state to the the 1-photon entangled NOON 
state shows that the introduction of entanglement does 
not remove the </> dependence of the Fisher information 
when arbitrary losses rx and Vy are present. However, 
when Vx = Vy" , the Fisher information, Fi(7r/2 — (/>), is 
independent of (p. This is in agreement with the result of 
Chen and Jiang, who derived the Fisher information for 
NOON state input using a master equation for a quantum 
continuous variable system for the case of symmetrical 
losses [6^. When losses are absent, — = 0, the 
Fisher information for input state given by Eq. ()64|) re- 
duces to F(0) = 1, independent of </>, as in the 1-photon 
Fock state without losses. 



2. 2-Photon NOON State Input into Lossy MZ 
Interferometer 



The 2-photon NOON state. 



|0) ^ [|20 ) + 102 )] 
(65) 
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rx 

FIG. 7. (Color) The Fisher information is plotted vs. for 
Tx ~ Ty, for Fock states (red) and NOON states (blue), for 
(/!> — 7r/4, where the Fisher information is non-zero for both 
Fock states and NOON states. 



has a peculiar behavior when input into a Mach-Zehnder 
interferometer with losses. The probabilities distribu- 
tions, given by Eqs. (l42t and (|63p . for the six measure- 
ment outcomes are independent of (j) and are given by 

P(20|(/., 20) = P(02|</., 20) = i{l-rl-rl + rlr^) 
P(11|0,2O) = 

P(1O|0, 20) - P(01|(/., 20) = ^{rl+rl- 2rlrl) 
Pi00\^,20) ^ rlrl 

(66) 



The measurement outcomes in Eq.(66) are independent 
of because of the Hilbert space geometry of the mea- 
surement operators, Il^(n,m), and input state vector in 
Eg. (1651) . see Eq. (j34p . For no prior information on the 
phase, p{(j)) — l/(27r), using Bayes' rule in Eq. (I?5]) . 
the phase probability densities are independent of (j), and 
are given by p(0|mn, 20) = l/(27r), for all measurement 
outcomes ^ = {m,n). Therefore, the 2-photon NOON 
state cannot be used in a Mach-Zehnder interferometer 
for determining the phase </>. However, such an arrange- 
ment can be useful in applications that require phase 
in-sensitive interferometry to be performed. In the limit 
of no loss, Tx — Ty = 0, the interferometer acts as a beam 
splitter and both photons come out the same port. 




FIG. 8. (Color) The Fisher information is plotted vs. for 
small dissipation in one arm, r-y — 0.10, for Fock states (red) 
and NOON states (blue), for = 7r/4. 

E. Imperfect Photon Number Detection 

Next, I consider the simplest example of imperfect 
photon-number detection. I assume that state prepara- 
tion is deterministic and that the interferometer is ideal, 
so there are no losses. I also assume that the input state 
is a pure state, in Eq. (ITT]) taking 

Psi^n = {', if 110) (67) 

^ ' [0, otherwise ^ ' 

where again I use the short- hand notation |10) for |1000). 
For this 1-photon input state, = |10), the no- loss 

Mach-Zehnder interferometer transfer matrix is given by 

( sin2<j!>, for j^^") = |10)and |V'°"') = |10) 
P7(V'°"*|V'"'» = I cos^fP, for jV''") = |10>and jV'""') = |01) 
0, otherwise 

(68) 

For the detection system, I assume that there is a prob- 
ability pd to detect the state correctly and a probability 
Px to detect the state incorrectly, where Pd + Px = ^- I 
am neglecting the possibility of an inconclusive measure- 
ment outcome. The matrix, Pd(^|-0°"*, 0), in Eq. (fTTj) 
describing the state detection is then given by 

pd(?i^°"*,0) = |^^; I^IZ (69) 

From Eq. (|lip . the probabilities P(^|0) for measurement 
outcomes are 

P {10\(j)) = Pd sin^ (t) + px cos^ (j) 
P{01\^)^pxsm^(f> + Pdcos^^ ^ > 
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FIG. 11. (Color) The fidelity (Shannon mutual information) 
is plotted as a function of dissipation rx , for — ry, for one- 
, two- and three-photon Fock state input, and for one- and 
two-photon NOON state input. 



FIG. 9. (Color) The Fisher information is plotted vs. 
for large dissipation in one arm, Vy = 0.90, for Fock states 
(red) and NOON states (blue), for (p — 7r/4. This example 
of high dissipation shows that the situation is complicated at 
high values of the dissipation in one arm and small values of 
dissipation in the other arm. 




-1 1 

(/> [radians] 

FIG. 10. (Color) The Fisher information is plotted vs. (j> for 
loss parameters = 0.3 and ry = 0.4, for NOON state input 
into a Mach-Zehnder interferometer for =1, 2, 3, 4, and 5 
photons. 



where ^ = (m, n) specifies that m and n photons are 
detected in output ports "c" and "d", respectively, see 
Eq. CiD-dini)- 



From Bayes' rule in Eq. 
probability density, 



(j25p . I find the conditional 
for the phase shift cj) for a 
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FIG. 12. (Color) The conditional phase probability density, 
p(<^|10), is plotted as a function of i\> for increasing detec- 
tor error probabilities -p^ = 0, 0.2, 0.4, 0.5, showing a loss of 
phase distinguishability. At pa; = 0.5, all phases are equally 
probable. 



given measurement outcome ^ to be 



(1 - -px) sin^ + cos^ I 
Pa; sin^(?!)-|- (1 -pa;)cos2i 



(71) 



Figure [T^] shows a plot of the phase probability density, 
p((/)|10) vs. 0, for different values of detector error proba- 
bility . With increasing probability of detecting the 
state incorrectly, the constrast in the phase probability 
decreases. Note that this contrast is not a "visibility" 
because p(0|lO) is a probability, and not an optical in- 
tensity. 

The fidelity (Shannon mutual information) defined in 
Eq. ^ is plotted in Fig. [T^l As expected for this sim- 
ple model, the fidelity H{M) decreases with increasing 
probability of incorrect detection, p^, reaching zero at 
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Px Detector Error Probability 

FIG. 13. (Color) The fidelity (Shannon mutual information) 
between the measurements and the phase is plotted vs. the 
probability of incorrect detection, p^. Consistent with the 
phase probability density plotted in Fig. 1121 the fidelity de- 
creases to zero at = 0.5 because there is no information on 
the phase in the measurements, so there is no discrimination 
between diflerent phases. 



Px — 0.5. Note that the fidelity is symmetric about 
Px = 0.5. 



IV. CONCLUSION 

I considered the experimentally relevant problem of 
determining the phase shift in one arm of a quantum 
interferometer when state creation is not perfectly deter- 
ministic, state propagation through the interferometer 
is non-unitary due to absorption losses in the interfer- 
ometer, and state detection is not ideal. In Section II, 
1 have argued that two types of information are useful 
for evaluating the quality of a parameter estimation de- 
vice, such as a quantum optical system used to determine 
phase shifts. First, fidelity (Shannon mutual informa- 
tion between measurements and parameter) is useful for 
deciding the overall quality of the optical system. The 
fidelity represents an average over probabilities of all pos- 
sible measurements and parameter values (phases). The 
fidelity is the metric to use when choosing or designing 
a system and the prior parameter (phase) distribution is 
unknown. Once a system is chosen, it is to be used in 
estimating the parameter based on measurements (data) , 
which is an estimation problem. At this point, the (clas- 
sical or quantum) Fisher information can be exploited. 



using the classical or quantum Cramer-Rao theorem, to 
estimate the variance of the parameter associated with 
its unbiased estimator. 

In Eq. pT|) . I have written down a general statisti- 
cal expression for the probability of a measurement out- 
come that simultaneously takes into account the three 
non-ideal aspects of real experiments: non-deterministic 
state preparation, losses in the interferometer, and non- 
ideal quantum state detection. This expression requires 
detailed models for each of the three non-ideal elements. 
In Section III, using simple, few-photon Fock states and 
NOON states, I give examples of applying Eq. ([TT|) . In 
subsection. A, of Section III, I look at a simple example 
of the effect of non-deterministic state creation, where 
there is a probability of creating one photon and a prob- 
ability of creating vacuum, as input into a Mach-Zehnder 
interferometer. As expected, the non-zero probability of 
creating a vacuum input leads to a probability for an in- 
conclusive measurement outcome, which in turn reduces 
the information on phase, as measured by Fisher infor- 
mation and fidelity (Shannon mutual information). 

In subsection B, of Section III, I have constructed a 
scattering matrix for a lossy (non-unitary) Mach-Zehnder 
interferometer. I find that for simple photon counting 
measurements, losses introduce a strong phase depen- 
dence in the classical Fisher information, making accu- 
racy of phase estimation dependent on the unknown true 
phase. In subsection, C and D, of Section III, I use the 
classical Fisher information and fidelity to examine in 
detail the propagation of Fock and NOON states, respec- 
tively, through the lossy Mach-Zehnder interferometer. 
Finally, in subsection E of Section III, I use Eq. (|lip . to 
look at the effect of imperfect photon number detection 
on determining the phase in a Mach-Zehnder interferom- 
eter with no losses, using the simplest model that takes 
into account a probability for incorrect detection of pho- 
ton number. 

The examples that 1 have used have been simple to 
illustrate the application of Eq. (fTTj). The theory can be 
applied to more complicated cases so that real experi- 
ments can be analyzed. 
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